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We summarize the analytical solution of the Chiral Perturbation Theory for the Hermitian 
Wilson Dirac operator of N c = 2 QCD with quarks in the fundamental representation. Results 
have been obtained for the quenched microscopic spectral density, the distribution of the 
chiralities over the real modes and the chiral condensate. The analytical results are compared 
with results from a Monte Carlo simulation of the corresponding Random Matrix Theory. 


1 Introduction 

The case of 5(7(2) QCD with fundamental quarks has attracted a great deal of interest in the 
past decades since it is a very interesting theory per se but it also shares many of the salient 
properties of SU (3) QCD, such as confinement and chiral symmetry breaking. An important 
difference between 577(2) and 577(3) is that 5(7(2) is pseudoreal, so that for an even number of 
flavors it is not hindered by the notorious sign problem at non zero chemical potential. Therefore, 
it is in an interesting playground for finite density studies, and one can study the phase diagram 
of 517(2) QCD in the T — /j plane 1,2 and get some insight in the qualitative properties of the 
phase diagram of QCD. One can get an understanding of the systematics but one needs to bear 
in mind that the two theories have quantitative differences. For 5(7(2) QCD the baryons are 
of bosonic nature and as a result this theory has a superfluid phase that is not encountered in 
ordinary QCD. Moreover, for beyond the Standard Model studies dedicated to the search of the 
conformal window, it has been shown that one can enter the conformal window with considerable 
less flavors compared to the 5(7(3) theory 3 . All these non-perturbative studies have been mainly 
through lattice simulations employing one of the most prominent discretizations, that of Wilson 
fermions. This discretization explicitly breaks chiral symmetry and consequently the low lying 
part of the Dirac spectrum gets affected in a very particular way. 

We study the discretization (cutoff) effects on the Dirac spectrum by employing the methods 
of Wilson Chiral Perturbation Theory (WchPT) and Wilson Random Matrix Theory (WRMT). 
This work extends previous work on 5(7(3) QCD 4 - 5 > 6 7,8,9 anc } we p 0 j n t Q ut the similarities and 
the differences among the two theories. 
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2 Wilson chPT 


Initially we concentrate on the study of the Hermitian Wilson Dirac operator D§ = 75 D\y where 
D\y is given by 

Dw = ^(V^ + Vp-|V M V*. (1) 

We work in the e— regime of WchPT where the mass m, the axial mass z and the lattice spacing 
a scale with the lattice volume such that the quantities 

mkS, zVS, and a 2 VW s (2) 


are kept fixed (S is the chiral condensate and Wg is the leading LEC). In this regime the 
partition function of WchPT factorizes and its mass dependence is given by a unitary matrix 
integral since the fluctuations of the zero momentum pions dominate completely the ones with 
non-zero momentum. The chiral partition function for Nf flavors reads 10 


ZN f (m, z ; a) = 


' SU (2Nf) / Sp(2N f ) 


d JJ e f SFTr((7+t/t)+f SVTr {U-U^)-a 2 VW 8 Tr (t/ 2 +t/t 2 ) ( 3 ) 


To access the spectral properties of the Dirac operator we employ the supersymmetric 
method. Adding an additional fermionic quark as well as an additional bosonic quark, with 
source terms in the mass as well as in the axial mass, one can differentiate with respect to these 
sources and compute the Green’s function of this theory in the standard way. By setting the 
source terms to zero one recovers the original partition function. In this proceeding we present 
our results graphically and we refer the reader to the original work for detailed derivations 9 . 

The first spectral observable that we consider is the quenched resolvent of defined as 




(4) 


whose discontinuity yields the spectral density of D§ 


P5(A) 


—IrnG^A + ie). 
7 r 


(5) 


The spectral density of D§ for N c = 2 exhibits a striking difference compared to N c = 3 
QCD. We see that there is a non-uniform convergence to the continuum limit 11 since there is 
a discontinuous jump by a factor of 1 /v /2 which is persistent for arbitrarily small values of the 
lattice spacing as it is clearly shown in Fig. 1. 

At the moment we have not been able to derive analytical results for the real eigenvalue 
density of D\y but we have obtained stringent upper and lower bounds of this distribution. The 
first distribution has been coined as the distribution of chiralities over the real modes of D\y 
and is defined as 

PxW = J2 5{\ + Ayv)sign[(Ayu|75|Arv}], (6) 

and together with 


PsCAs = 0 ,m;a) 



/ 6(\w + m) \ 

\ k*m*>i r 


(7) 


we have the following inequality for the density p re ai(A) of the real modes of D\y, 


Px( A) < /Oreal(A) < /0 5 ( A 5 = 0,m = A). (8) 

In Fig. 2 we see that for large values of the lattice spacing these distributions are bounding the 
real eigenvalue density from above and below. While for values of a close to zero we see that p x 
is a very good approximation of p rea i . 





Figure 1 - The spectral density of D 5 for z = 0, m = 0 and values of a as in the legend of the figures. We see that 
the spectral density at zero virtuality jumps by a factor l/y/2 when 0. The solid curves denote the analytical 
results while the histograms represent the data of a simulation of 10 s 200 x 200 random matrices. On the RHS 
we see the discontinuous behavior of ps for very small values of a. 




Figure 2 - In the left figure we show the spectral density of D$ for v = 1, m = 2 and a = 0.025, and in the right 
figure we display the spectral density of D 5 for v = 1, m = 2 and a = 0 (black), a = 0.15 (green) and a = 0.3 
(red). The numerical results (histogram) are compared to the analytical result (curve). Also shown are the results 
for a = 0 (black curve) with a <5-peak at A 5 EI/ = 2. The S peak corresponding to the zero modes gets broadened 
with a width proportional to the lattice spacing a. It is interesting to see that for larger values than a = 0.3 the 
peak of the zero modes has almost completely dissolved. 




Figure 3 - The distribution of the chiralities (blue curves) and the inverse chiralities (red curves) over the spectrum 
of D 5 compared to Monte Carlo simulations for v = 1, a — l(left) and a = 0.125(right). The density p Tea i of 
the real eigenvalues of D\y (black histogram) is always bounded by p x and ps(As = 0) and can be accurately 
approximated by p x for small values of the lattice spacing. The analytical results are in agreement with Monte- 
Carlo simulations (blue and red histograms). 
























3 Conclusions 


We have studied the microscopic spectrum of the Wilson Dirac operator for QCD with two colors 
and quarks in the fundamental representation. The discretization effects for SU( 2) QCD are 
mostly similar to the case of SU( 3) QCD. As we have seen from the figures, the main effects are 
the broadening of the topological peak and the tail states that enter the gap of the Hermitian 
Dirac operator. Also in this case the gap closes when entering the Aoki phase. Similarly to the 
case of ordinary QCD additional real eigenvalues could potentially be an issue in the context of 
SU (2) QCD. These originate from the collision of complex conjugate pairs which eventually enter 
the real axis, and on average, we have a u+1 additional real eigenvalues when we are restricted 
sufficiently close to the continuum limit. However, note that the suppression of these modes 
(for v ^ 0) is not as strong as for ordinary QCD 6 . Another noteworthy property also found for 
QCD with three colors is that based on the Hermiticity properties of the Dirac operator one can 
pinpoint 9 the sign of Wg. In our conventions when taking into account only Wg in our analysis 
we conclude that it should have a positive sign. The sign of Wg is important because it controls 
if the theory enters or not the Aoki phase. Our results provide an analytical handle on the 
smallest eigenvalues of the Wilson Dirac operator which could potentially seriously compromise 
the Monte Carlo simulation if they become almost zero. One needs to point out that the absence 
of level repulsion from the origin for the case of two-color QCD makes this effect much more 
severe compared to 577(3) QCD. 
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